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I. INTRODUCTION 


A previous  technical  report  [lj  discussed  the  use  of  a logistic 
function  in  development  of  an  impact  acceleration  injury  prediction  model 
based  on  empirical  data.  Another  report  [2]  described  a Monte  Carlo  study 
conducted  to  assess  the  accuracy  with  which  model  parameters  and  injury 
probabilities  could  be  estimated. 

The  model  under  consideration  is  based  on  the  assumption  of  an 
underlying  functional  relationship  of  the  form 


P(*) 


{l  + exp[-($Q 


k 

+ ze4x  )]}' 
1 1 1 


(i) 


where : 

x * (x^,...,x,)  denotes  the  set  of  independent  variables  considered, 

denotes  a set  of  parameter  values, 

and  P(x)  denotes  the  true  probability  of  Injury  corresponding  to  x . 

This  report  considers  the  application  of  such  a model  to  observed 
data  from  a set  of  twenty-eight  -G^  accelerator  runs  Involving  subhuman 
primates  (Rhesus  monkeys)  with  securely  restrained  torso  and  unrestrained 
head.  The  data  was  collected  by  the  Naval  Aerospace  Medical  Research 
Laboratory  (NAMRL)  Detachment  as  part  of  its  research  effort  on  accelera- 
tion impact  Injury  prevention. 

Two  prediction  models  were  constructed  from  this  data,  one  based  on 
head  dynamic  response  only  and  the  other  based  on  sled  acceleration  profile 
only.  The  first  model  was  derived  from  three  variables  distilled  from 
head  dynamic  response  time  trace  data.  These  three  variables  were: 
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(1)  peak  head  angular  acceleration  (resultant)  measured  In  radians/ 

2 

sec  , 

(2)  peak  head  linear  acceleration  (resultant)  measured  In  meters/ 

2 

sec  , 

and  (3)  peak  head  angular  velocity  (resultant)  measured  In  radians/sec  . 
The  second  model  was  based  on  two  variables  describing  sled  acceleration: 

(1)  peak  sled  acceleration  measured  In  G's 
and  (2)  rate  of  sled  acceleration  onset  measured  in  G/sec. 

Surprisingly,  the  model  based  on  sled  profile  variables  provided 
better  predictions  than  the  model  based  on  head  dynamic  response  variables. 
Because  of  this,  both  sets  of  variables  were  combined  into  one  overall 
five-variable  set,  which  was  then  used  In  development  of  a prediction  model. 
It  was  found  that,  for  the  relatively  small  amount  of  data  available, 
inclusion  of  peak  sled  acceleration  and  any  one  of  the  other  variables 
resulted  in  a model  which  yielded  predictions  in  almost  perfect  agreement 
with  the  observations. 
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Figure  1:  The  Data  Set 


may  prove  unimportant  and  should  therefore  not  be  Included  In  a final 
model.  To  judge  the  contribution  of  variables,  likelihood-ratio  tests 
may  be  conducted.  This  procedure  may  be  used  in  conjunction  with  "nested" 
models. 

In  the  present  context,  one  model  will  be  said  to  be  nested  within 
another  if  the  second  model  contains  all  variables  of  the  first  model 
plus  one  or  more  additional  variables.  Thus,  a model  which  contained 
variables  x^  , and  x^  would  be  nested  within  a model  which  contained 
only  x and  x . 


1 2 

To  test  a hypothesis  that  a model  containing  variables  (x^ x^^) 

is  a significant  improvement  over  a model  containing  variables  (x^,... 
a Chi-square  statistic  may  be  used.  The  procedure  is  to  calculate: 

“ -2  log  likelihood  for  model  containing  (x^, . . . ,x^) 
and  L - -2  log  likelihood  for  model  containing  (x  , ...,x ). 

2 1 iCTtn 

Under  the  null  hypothesis  that  the  m additional  variables  ) 

do  not  result  in  an  improved  model,  the  statistic  has  an  approximate 

Chi-square  distribution  with  m degrees  of  freedom.  Thus,  the  hypothesis 
may  be  tested  by  comparing  the  value  of  - L2  with  the  upper  percentage 
points  of  the  appropriate  Chi-square  distribution. 


B.  HEAD  DYNAMIC  RESPONSE  VARIABLES 

As  previously  mentioned,  three  head  dynamic  response  variables  were 
considered  in  model  development.  In  the  following  discussion  these  vari- 
ables will  be  denoted  by  Xj,  x^,  and  x^,  where: 

x^  is  peak  head  angular  acceleration  (reaultant)  measured  in  radians/ 
2 

sec  , 
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x is  peak  head  linear  acceleration  (resultant)  measured  in  meters/ 

2 

sec^  , 

and  x^  is  peak  head  angular  velocity  (resultant)  measured  in  radians/sec  . 
The  first  stage  in  model  construction  based  on  variables  x^,  x^  and  x^ 
involved  examination  of  all  possible  models  including  these  variables  or 
any  subset. 

For  a given  number  of  variables,  that  model  which  yielded  the  smallest 
-2  log  livelihood  value  was  se1 ected  as  the  ''best".  As  can  be  seen  from 
Figure  2,  the  best  one-variable,  two-variable,  and  three-variable  models 
are  those  based  on,  respectively, 

(1) 


(2)  *2.  x3 


(3)  x , x 2»  x 3 • 


Because  of  the  nesting  in  these  models,  the  relative  contribution  of 

variables x , x , and  x may  be  tested  in  that  order. 

2 3 1 

Figure  3 presents  a summary  of  the  relevant  test  procedure.  In 

the  first  stage,  x^  was  tested  to  determine  whether  it  significantly 

improved  a model  which  assumed  constant  probability  over  all  values  of  the 

three  head  dynamic  response  variables.  The  observed  Chi-square  value  of 

11.26,  which  is  statistically  significant  at  the  .0008  level,  indicated 

that  this  variable  did  result  in  an  improved  model. 

The  second  stage  of  testing  involved  consideration  of  the  addition 

of  another  variable  to  the  model  which  included  only  variable  x^  . 

Because  the  best  two-variable  model  was  based  on  x and  x . the  effect 

2 3 
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Variable  Set 

-2  Log  Likelihood 

Constant  Only 

29.10 

X 

19.30 

1 

X2 

17.84 

X3 

19.89 

. 

V x2 

17.82 

Xl*  x3 

19.07 

*2»  x3 

17.47 

V V X3 

16.91 

denotes  peak  head  angular  acceleration 
denotes  peak  head  linear  acceleration 
x^  denotes  peak  head  angular  velocity 

Figure  2:  Head  Dynamic  Response  Variable  Sets  and 
Associated  -2  Log  Likelihood  Values 


Test  1:  X£  against  Constant  Only 

L - L * 29.10  - 17.84  - 11.26  (1  d.f.,  p = .0008) 

Significant 


Test  2:  (x2,  x^)  against  X2 

L - L2  - 17.84  - 17.47  - 0.37  (1  d.f.,  p > .50) 

Nons ignif leant 

Test  3:  (x  , x2,  x3)  against  x2 

L - L - 17.84  - 16.91  = 0.93  (2  d.f.,  p > .50) 

Nonsignificant 

x^  denotes  peak  head  angular  acceleration 

x denotes  peak  head  linear  acceleration 
2 

denotes  peak  head  angular  velocity 


Figure  3:  Testing  the  Significance  of  the  Head  Dynamic 
Response  Variables 
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m 

f 

tm 

If 

of  including  the  latter  variable  was  examined.  The  addition  of  x^  to  the 

model  resulted  in  an  observed  Chi-square  value  of  0.37,  which  is  not 

statistically  significant.  Likewise,  a model  based  on  all  three  variables 

(x^,  x^,  x^) , when  tested  against  the  model  based  on  X2  only,  resulted  in 

a statistically  nonsignificant  Chi-square  value  of  0.93. 

Thus,  based  on  the  set  of  data  under  consideration,  a predictive 

model  which  includes  only  variable  x^  appears  to  provide  the  best  results. 

It  is  interesting  to  note  that  x^  is  a linear  component,  while  both  x^ 

and  x^  are  angular  components.  This  provides  some  tentative  empirical 

evidence  that,  for  -G  acceleration,  the  primary  component  correlating 

x 

with  Impact  acceleration  injury  (more  correctly,  fatality)  may  be  linear 
rather  than  angular.  However,  this  evidence  is  far  from  concl  ive 
because  of  the  small  data  base  and  the  high  degree  of  intercorrelations 
between  the  observed  values  of  the  variables. 

The  resulting  model  is  given  by: 

A -1 

P(x2)  - (1  + exp [-(-4. 344  + 0.001369  x^]}  (2) 

A 

where,  of  course,  P(x^)  denotes  the  fitted  model.  Figure  4 presents,  for 
this  model,  a comparison  of  observed  probability  (i.e.,  0 or  1)  and  pre- 
dicted probability,  where  the  observations  are  arranged  in  order  of 
Increasing  predicted  probability.  As  can  be  seen  from  this  figure,  there 
is  reasonable  agreement  between  predictions  and  observations. 

C.  SLED  PROFILE  VARIABLES 

A procedure  similar  to  that  described  in  the  previous  section  was 
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0 
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LX  1899 

A04101 

0 
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580.0 

LX  189  0 

A04101 

0 

0.0294 

618.0 

LX1084 

A0  392 1 

0 

0.0298 

630.0 

LX  1085 

A03921 

0 

0.0300 

635.0 

LX1086 

A03921 
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0.0302 
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A03921 
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0.0304 
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A03921 
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A04101 

0 
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A03921 
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A03921 

0 
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A04101 
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1540.0 
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A03948 
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0.  1092 
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A04101 
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0.1153 

1685.0 

LX1902 

A04101 

0 

0. 1285 
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LX  1903 

A04101 
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0.1316 

1795.0 

LX1363 

AC3935 

1 

0. 1569 

1945.0 

LX  1891 

A03943 

0 

0.1587 

1955.0 

LX1893 

AC3924 

0 

0.2576 

2400.0 

LX  18  96 

AO  3946 

1 

0.2805 

2485.0 

LX  18 94 

A03933 

0 

0. 2959 

2540.0 

LX  1359 

A04099 

0 

0.3686 

2780.0 

LX1905 

A04101 

1 

0. 3814 

2820.0 

LX1895 

A03951 

1 

0.3912 

2850.0 

LX1362 

A03935 

0 

0.7806 

4100.0 

LX  1360 

A04099 

1 

0.9421 

5210.0 

LX1365 

A0  3921 
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0.9998 

9250.0 

Figure  4:  A Comparison  of  Observed  and  Predicted 
Probability  for  Model  Based  on  Peak 
Head  Linear  Acceleration 
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applied  to  two  sled  profile  variables  which  in  the  following  discussion 


are  denoted  by  z 


1 


and  z^. 


where 


z 

1 

and  z 

2 


is  peak  sled  acceleration  measured  In  G's 

is  rate  of  sled  acceleration  onset  measured  in  G/sec  . 


As  Figure  5 Indicates,  the  "best"  one-variable  model  is  based  on  z^  . 

Figure  6 provides  a comparison  of  predicted  and  observed  probability 
for  this  model,  which  is  given  by 

P(zx)  - U + exp [-(-49. 81  + 0.4472  z^]}"1  . (3) 


As  can  be  seen  from  this  figure,  the  agreement  between  predictions  and 
observations,  although  not  perfect,  is  relatively  good. 

Figure  7 summarizes  the  testing  of  the  model  improvement  resulting 
from  the  inclusion  of  variables  z^  and  z^  in  that  order.  The  results  of 
the  test  in  the  first  stage  (observed  Chi-square  value  of  24.34,  statis- 
tically significant  at  the  .0001  level)  indicated  that  variable  z^  provided 
an  improved  model.  Likewise,  the  second  stage  of  testing  showed  that  Z2 
should  be  Included  in  the  model  also,  since  the  observed  Chi-square  value 
was  4.75,  which  is  statistically  significant  at  the  .029  level. 

The  resulting  fitted  model  based  on  both  variables  z^  and  z^  is 
given  by 

A -1 

P(Zl,z2)  - (1  + exp [-(-4463.0  + 38.88  zx  + .01816  z2>j)  (4) 

A comparison  of  observed  0/1  probabilities  with  predicted  probabilities 
la  given  in  Figure  8.  As  indicated  in  this  figure,  the  agreement  between 
observed  and  predicted  probability  is  almost  perfect. 
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Variable  Set  -2  Log  Likelihood 


Constant  Only  29.10 


4.76 

13.92 


z^  denotes  peak  sled  acceleration 

denotes  rate  of  sled  acceleration  onset 


Figure  5:  Sled  Acceleration  Profile  Variable  Sets 
and  Associated  -2  Log  Likelihood  Values 
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Figure  6:  A Comparison  of  Observed  and  Predicted 
Probability  for  Model  Based  on  Peak 
Sled  Acceleration 
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Test  1:  z against  Constant  Only 

L - L - 29.10  - 4.76  - 24.34  (1  d.f.,  p = .0001) 

Significant 


Test  2:  (z^,  Z2>  against  z^ 

L - L - 4.76  - 0.01  - 4.75  (1  d.f.,  p - .0293) 

Significant 


z^  denotes  peak  sled  acceleration 


z^  denotes  rate  of  sled  acceleration  onset 


Figure  7: 


Testing  the  Significance  of  the  Sled 
Acceleration  Profile  Variables 
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Figure  8:  A Comparison  of  Observed  and  Predicted 
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D.  COMBINED  HEAD  AND  SLED  VARIABLES 


In  the  two  previous  sections  the  three  head  dynamic  response  variables 
and  the  two  sled  profile  variables  were  regarded  as  defining  two  separate, 
independent  sets  of  variables.  In  this  section,  these  five  variables  are 
considered  as  comprising  one  overall  set  from  which  prediction  models  are 
developed . 

As  Figure  9 indicates,  z^  (peak  sled  acceleration)  provided  the  best- 
fitting one-variable  model.  From  this  figure  it  can  also  be  seen  that 
four  models  are,  for  all  practical  purposes,  tied  in  the  competition  for 
the  best  two-variable  model.  These  are  the  models  based  on  (zp  x^), 

(z  , x2),  (zp  x3),  and  (Zp  z2> . 

In  general,  the  testing  procedure  involves  a test  of  z^  against  a 
constant  probability  model  followed  by  a test  of  whether  the  addition  of 
a second  variable  provides  significant  improvement.  Because  the  -2  log 
likelihood  values  for  the  addition  of  any  one  of  the  four  variables  Xp 
x2>  x^,  or  z2  to  the  model  are  essentially  the  same  (0.00  or  0.01),  the 
test  results  are  equivalent  to  those  illustrated  for  z2  in  Figure  7. 

By  examining  Figure  8 and  Figures  10  through  12,  it  can  be  seen  that 
any  of  the  four  variable  pairs  result  in  a model  which  yields  predictions 
in  almost  perfect  agreement  with  the  observed  data.  The  model  involving 
z^  and  z2  is  given  by  (4)  in  the  previous  section.  The  three  remaining 
models  are: 

~ -1 

P(zr  xx)  - {1  + exp [-(-354.0  + 2.767z1  + 0.00108x^1)  (5) 

P(zp  x2)  - {1  + exp [-(-196. 7 + 1.622Zl  + 0.00345x2)] )_1  (6) 

P(zp  x3)  - (1  + exp[- (-248.6  + 2.009zi  + 0.11938x3) J }_1  (7) 
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x^  denotes  peek  heed  enguler  velocity 
z]T  denotes  peek  sled  ecceleretion 
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Figure  9:  Heed  end  Sled  Acceleretlon  Verleble  Sets 
end  Assocleted  -2  Log  Likelihood  Velues 
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Figure  12:  A Comparison  of  Observed  and  Predicted 
Probability  of  Model  Baaed  on  Peak  Sled 
Acceleration  and  Peak  Head  Angular  Velocity 


Velocity 


In  summary.  It  may  be  stated  that,  based  on  the  data  set  considered. 


i 


I 


I 


peak  sled  acceleration  is  the  single  most  useful  variable  In  making 
accurate  predictions.  However,  its  effectiveness  can  be  significantly 
improved  by  adding  any  one  of  the  four  other  variables: 

(a)  z^,  rate  of  sled  acceleration  onset 

(b)  x^,  peak  head  angular  acceleration 

(c)  x , peak  head  linear  acceleration 

2 

(d)  x^,  peak  head  angular  velocity. 

E.  PREDICTION  OF  CRITICAL  ENVELOPES 

From  a fitted  model  such  as  (2)  through  (7),  a critical  envelope  can 
be  predicted.  This  envelope  defines  those  combinations  of  independent 
variables  for  which  the  predicted  probability  of  injury  (or  fatality)  is 
greater  than  some  given  amount.  In  the  present  situation,  suppose  it 
were  desired  to  restrict  the  variable  values  to  a region  in  which  the 
probability  of  fatality  were  less  than  some  small  probability  Pq  (such 

A 

as  .01  or  .05).  In  other  words,  the  predicted  probability  P(x)  would  be 

less  than  P_  . From  this  it  follows  that: 

0 

A 

P(x)  < PQ 

A ICa  _ i 

{1  + exp[-(eQ  + EB^)]}  < ?0 

* k* 

B0  + E&iXi  < ln[P0/(l  - PQ)J  . 

From  the  fitted  head  dynamic  response  model  given  by  (2),  the 
critical  envelope  at  P^  ■ .05  is 
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-4.344  + .001369x2  < -2.944 
or  equivalently, 
x2  < 1022.6  . 

Thus,  the  probability  of  fatality  is  predicted  to  be  less  than  5%  if 
peak  head  linear  acceleration  is  less  than  1022.6  radians/sec  . This  is 
shown  graphically  in  Figure  13. 

Similarly,  the  critical  envelope  at  - .05  based  on  sled  profile 
variables  incorporated  into  fitted  model  (4)  is 

-4463.0  + 38.80z1  + .01816z2  < -2.944 

or  equivalently, 

2136. 56Zi  + z2  < 245,598.01  . 

This  is  shown  graphically  in  Figure  14.  The  5Z  critical  envelopes  for 
the  other  three  variable  pairs  are  illustrated  in  Figures  15  through  17. 
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Figure  13:  Boundary  for  Predicted  Fatality  Probability  of  5X 
(x  Denotes  Fatality,  0 Denotes  Nonfatality) 


III.  SUMMARY  AMD  DISCUSSION 


Because  of  differing  Initial  head  position  of  the  experimental 
subjects.  It  was  postulated  a priori  that  the  sled  acceleration  profile 
would  yield  less  sensitive  Independent  variables  than  head  dynamic 
response  would.  Using  a common  data  base,  two  different  models  were 
constructed,  one  based  on  sled  profile  variables  and  the  other  based 
on  head  dynamic  response  variables.  Although  the  latter  model  provided 
a reasonable  fit  given  the  small  size  of  the  data  set,  the  other  model 
(based  on  sled  profile  variables)  resulted  In  a much  better  fit.  This 
can  be  seen  by  comparing  Figures  4 and  8. 

Since  it  is  Intuitive  that  a model  based  on  head  dynamic  response 
should  provide  predictions  which  are  as  good  as  or  better  than  those 
from  a model  based  on  sled  profile,  some  explanation  is  required.  There 
are  a few  possible  reasons  for  this  anomalous  result.  It  may  be  that 
the  wrong  variables  were  extracted  from  the  head  dynamic  response  time 
traces,  and  other  variables  would  have  more  correctly  conveyed  the  Infor- 
mation within  these  time  traces.  On  the  other  hand,  the  correct  variables 
may  have  been  selected  but  errors  may  have  been  present  in  their  measure- 
ment. In  addition.  It  is  possible  that  the  small  sample  size  resulted  In 
a spurious  result. 

It  must  be  realized,  however,  that  consideration  of  any  of  the  three 
head  dynamic  response  variables  In  conjunction  with  peak  sled  acceleration 
provides  a perfectly  fitting  model.  Thus,  there  la  evidence  that  head 
dynamic  response  variables  and  sled  profile  variables  may  be  used  together 
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to  provide  good  results.  Nonetheless,  It  Is  still  an  open  question 
of  why  head  dynamic  response  variables  alone  did  not  provide  a model 
which  performed  as  well  as  the  one  based  on  sled  profile  variables  alone. 

In  any  event,  additional  accelerator  runs  are  warranted,  particularly 
In  the  region  defined  by 


100  < zx  < 125 


10,000  < z2  < 20,000 


where,  as  previously  defined. 


z^  Is  peak  sled  acceleration  measured  in  G's 


and  z^  is  rate  of  sled  acceleration  onset  measured  In  G/sec  . 


This  should  result  In  valuable  information  around  the  apparent  boundary 
between  fatality  and  nonfatality. 
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